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ABSTRACT 

Using published measurements of dielectronic recombination (DR) resonance 
strengths and energies for C V to C IV and O VIII to O VII, we have calculated 
the DR rate coefficient for these ions. Our derived rates are in good agreement with 
multiconfiguration, intermediate-coupling and multiconfiguration, fully-relativistic 
calculations as well as with most LS coupling calculations. Our results are not in 
agreement with the recommended DR rates commonly used for modeling cosmic 
plasmas. We have used theoretical radiative recombination (RR) rates in conjunction 
with our derived DR rates to produce a total recombination rate for comparison with 
unified RR+DR calculations in LS coupling. Our results are not in agreement with 
undamped, unified calculations for C V but are in reasonable agreement with damped, 
unified calculations for O VIII. For C V, the Burgess general formula (GF) yields a rate 
which is in very poor agreement with our derived rate. The Burgess & Tworkowski 
modification of the GF yields a rate which is also in poor agreement. The Merts et 
al. modification of the GF yields a rate which is in fair agreement. For O VIII the 
GF yields a rate which is in fair agreement with our derived rate. The Burgess &: 
Tworkowski modification of the GF yields a rate which is in good agreement. And 
the Merts et al. modification yields a rate which is in very poor agreement. These 
results suggest that for AA'^ = 1 DR it is not possible to know a priori which formula 
will yield a rate closer to the true DR rate. We describe the technique used to obtain 
DR rate coefficients from laboratory measurements of DR resonance strengths and 
energies. For use in plasma modeling, we also present easy-to-use fitting formulae for 
the experimentally derived DR rates. 



Subject headings: atomic data — atomic processes 



1. Introduction 



Carbon and oxygen are two of the most abundant elements in the universe and lines from 
these elements provide valuable plasma diagnostics for almost all classes of cosmic sources. 
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Essential for many of these diagnostics are accurate electron-ion recombination rate coefficients, 
particularly of dielectronic recombination (DR), which for most ions in electron-ionized plasmas 
is the dominant means of electron capture (Arnaud & Rothenflug 1985). Producing accurate 
theoretical DR rate coefficients is, however, theoretically and computationally challenging. 

In the past, semi-empirical expressions such as the Burgess (1965) formula along with 
modified versions by Burgess & Tworkowski (1976) and Merts et al. (1976) were developed 
to calculate DR rates. More recently, a number of more sophisticated theoretical approaches 
have been used to calculate DR, among them single-configuration LS'-coupling (Bellantone Sz 
Hahn 1989), multiconfiguration intermediate-coupling (Pindzola, Badnell, &; Griffin 1990), and 
multiconfiguration fully-rclativistic (Chen 1988) techniques, as well as undamped and damped, 
unified radiative recombination (RR) and DR calculations in L^-coupling (Nahar & Pradhan 
1997; Nahar 1999). Approximations, though, need to be made to make any of these techniques 
computationally tractable (Hahn 1993). Currently, sophisticated DR calculations are non-existent 
for many ions, and in the absence of anything better, semi-empirical formulae are often still used 
for plasma modeling. 

Laboratory measurements can be used to test the different theoretical and computational 
techniques for calculating DR. Recently, Savin et al. (1997, 1999) developed a technique for 
obtaining rate coefficients from laboratory measurements of DR resonance strengths and energies. 
They successfully used this technique to derive rates for An = DR of Fe XVIII and Fe XIX and 
to benchmark existing theoretical calculations. Here, we describe this technique in detail for the 
first time and apply it to recent DR measurements in C V and O VIII. Kilgus et al. (1990, 1993) 
and Mannervik et al. (1997) have measured the resonance strengths and energies for DR of C V 
to C IV and O VIII to O VII. We use their results to produce DR rate coefficients to benchmark 
existing C V and O VIII DR calculations and to provide rates for use in plasma modeling. 

In electron-ionized plasmas, lines from heliumlike C V and hydrogenic O VIII trace gas at 
Te ~ 1051-5-9 K and ~ lO*'-^ K, respectively (Arnaud & Rothenflug 1985; Mazzotta et al. 1998). 
C V and O VIII lines have been observed in solar spectra (Doschek & Cowan 1984) and O VIII 
lines in supernova remnants (Winkler et al. 1981). And with the upcoming launches of Chandra 
and XMM and the high-resolution spectrometers aboard, C V and O VIII lines are expected to be 
seen in may other electron-ionized, cosmic sources. 

Using different heavy- ion storage rings, Kilgus et al. (1993) and Mannervik et al. (1997) have 
measured DR for C V via the capture channels 

C4+(is2) + e- ^ C''+{ls2lnl') (n = 2, . . . , n„,,). (1) 

where nmax ~ 28 for the results of Kilgus et al. and 16 for the results of Mannervik et al. Kilgus 
et al. (1990) have also measured DR for O VIII via the capture channels 

0^+(ls) + e- ^ 0^+{2lnl') (n = 2, . . . , nmax) (2) 

where nmax ~ 69. The radiative stabilization of these autoionizing C V and O VII states to 
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bound configurations results in DR. Details of the experimental techniques used are given in the 
references cited. 

The paper is organized as follows: We describe in Section |2| how one produces a DR rate 
coefficient using measured DR resonance strengths and energies. In Section ^ we present the 
resulting rate coefficients and compare the derived DR rates with published theoretical rates. We 
also give a simple fitting formula for use in plasma modeling. 



2. Method of Calculation 

DR is a resonance process consisting, in the zero-density limit, of an infinite number of 
resonances. The DR rate coefficient a for a plasma with a Maxwell-Boltzmann electron distribution 
is given by 

a{Te)= [Y.ad{E)ve{E)P{Te,E)dE (3) 

where Tg is the electron temperature; (Td{E) is the energy-dependent DR cross section for resonance 
d; Ve{E) is the relative electron- ion velocity at energy E, which is taken to be the electron energy 
as the ions are nearly to stationary in the center-of-mass frame; and the sum is over all DR 
resonances. The Maxwell-Boltzmann distribution P{Te,E) is given by 

2/71/2 / _p \ 

where ks is the Boltzmann constant. 

Kilgus et al. (1990, 1993) and Mannervik et al. (1997) published measured DR resonance 
strengths ad and energies E^. The DR resonance strength is defined 

rE+5E/2 

CTd{E) = / ad{E')dE' (5) 

JE-5E/2 

where ad is the cross section for a resonance or group of resonances labeled d and {E — 5E,E + 6E} 
is a region in energy chosen such that it contains only those resonances comprising d. 

Here we are interested in calculating rate coefficients. This involves convolving the DR 
resonances with the slowly varying function P(Te,E). Because the energy widths of the measured 
resonances are smaller than the scale over which P(Te,E) changes, for our purposes we can 
accurately approximate ad{E) as 

ad{E) = ad{Ed)6{E - Ed) (6) 

where Ed is the energy of resonance d and 6{E — Ed) is the Dirac delta function. The DR rate 
coefficient for Maxwellian plasmas is found by substituting Equation ^ into Equation ^ which 
yields 

a{Te) = ^ dd{Ed)ve{Ed)P{T„Ed). (7) 
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Kilgus et al. (1993) and Mannervik et al. (1997) do not report measured resonance energies 
for capture by C V into levels where n > 4. To calculate these resonance energies En we use the 
Rydberg formula 

En = AE- q^Ro^/n^ (8) 

where g = 4 is the charge of the ion before recombination, IS.E = 307.8 eV is the energy of the 
Is^Q-Sq) — \s2pQ-Pi) core excitation (Kelly 1987), and R^o is the Rydberg constant. For O VIII, 
q = 7 and AE = 653.6 eV is the energy of the 15(^5*1/2) — 2^(^^3/2,1/2) core excitation (Kelly 
1987). 

Mannervik et al. (1997) estimate that they measured DR for capture into levels rimax ^ 16. 
Kilgus et al. (1993) do not report a value of Umax, so we derive an estimate here. Using their ion 
energy, the bending radius of the dipole magnets in their experiment (Linkemann 1995), and the 
semiclassical formula for field ionization (Brouillard 1983), we estimate electrons captured into 
levels where n = ricut ^ 19 will be field ionized and thus not detected. However, the captured 
electrons can radiatively decay below ricut during the ~ 5.1 m distance they travel between the 
electron- ion interaction region in the experiment and the dipole magnet (cf., Kilgus et al. 1992; 
Savin et al. 1997). Using the ion velocity, the hydrogenic formula for radiative lifetimes of Marxer 
& Spruch (1991), and calculations which show that DR for heliumlike ions essentially populates 
only angular momentum levels where / < 4 (Chen 1986), we estimate that electrons captured into 
levels Umax ^ 28 will radiatively decay below ricut and thus were detected by Kilgus et al. For 
OVIII, Kilgus et al. (1990) estimate Umax ~ 69. 

For a given DR series, as n increases the energy spacing between DR resonances decreases. 
Due to the energy resolutions of the experiments, above some n level it is not possible to resolve 
the DR resonances and determine individual values of an- For these high n levels, Kilgus et al. 
(1990; 1993) presents total resonance strengths summed from n > 9 to rimax for C V and from 
71 > 8 to Umax for O VIII. We divide and spread out these summed resonance strengths into ~ 1 
eV wide bins when using Equation |7[ The "resonance" energies of these bins are chosen to lie 
between Eg and En^^^ for C V and between and En^^^ for O VIII. The bin widths are smaller 
than the scale over which P{Te, E) changes and the small errors in the exact resonance energies of 
these high n levels has an insignificant affect on the derived DR rates. 

3. Results and Discussion 

In Figures [l| and ^ we present the C V and O VIII DR rate coefficients, respectively, derived 
using the measured resonance strengths and energies of Kilgus et al. (1990, 1993) and Mannervik 
(1997). The unmeasured contribution to the DR rate due to capture into levels where n ^ Umax 
is predicted to be insignificant (Chen 1986; Pindzola et al. 1990). The absolute uncertainties in 
the derived DR rates are estimated to be ^ ±25% which corresponds to the reported absolute 
experimental uncertainties. 
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Existing theoretical DR rates are also shown in Figures |T| and |2| as well as the RR rates of 
Verner & Ferland (1996). For C V, the single-configuration L5-coupling calculations of Bellantone 
& Hahn (1986) and Romanik (1988), multiconfiguration intermediate-coupling calculations of 
Badnell, Pindzola, & Griffin (1990), and multiconfiguration fully-relativistic calculations of 
Chen (1988) are all in good agreement with our experimentally inferred rate. The L5-coupling 
calculations of Younger (1983) are ~ 30% larger than the experimental rate. The recommended 
and commonly used DR rate of Shull & van Steenberg (1982) peaks at a value ~ 45% larger than 
ours and has a steeper low Tg behavior. 

Pindzola et al. (1990) calculated resonance strengths for O VIII using intermediate-coupling, 
multiconfiguration Hartree-Fock (MCHF) and intermediate-coupling, multiconfiguration Thomas- 
Fermi (MCTF) techniques. Using their MCHF results for capture into levels n < 6 and their 
MCTF results for n > 7, as well as the experimental resonance energies, we calculate the 
corresponding DR rate using Equation |^. The resulting rate is in good agreement with our derived 
rate. The single-configuration LS'-coupling rate of Bellantone & Hahn (1986) is also in good 
agreement with the experimental rate, though with decreasing Tg their rate falls off sooner than 
ours. The recommended and commonly used DR rate of Shull & van Steenberg (1982) is ~ 31% 
larger than our derived rate. Not shown in Figure |^, for reasons of clarity, is the rate of Zhdanov 
(1978) whose peak rate is ~ 7 times larger than the experimental rate. 

Nahar & Pradhan (1997) and Nahar (1999) present unified RR-I-DR rates. To compare with 
their results we add the RR rates of Verner & Ferland (1996) to the experimentally derived DR 
rates. Thus we treat RR and DR as independent processes and do not allow for the possibility of 
interference between the two recombination channels. The validity of this approach is supported 
by recent theoretical work (for DR on systems ranging in complexity from C II and Mg II to U 
XCII) which has shown the effect of interferences to be small (Pindzola et al. 1992). Figures H 



and ^ show for C V and O VIII, respectively, the unified RR-I-DR rate of Nahar & Pradhan and 
the total RR-I-DR rate using the derived DR results and the RR results of Verner & Ferland. For 
C V, at peak value the undamped, unified rate of Nahar &: Pradhan is ~ 43% larger than our 
resulting total RR-I-DR rate. This is consistent with the estimate by Pradhan & Zhang (1997) 
that allowance for radiation damping would reduce the undamped rate by ~ 20 — 30%. For O 
VIII, at peak value the damped, unified rate of Nahar is ~ 20% larger than our resulting total 
RR-hDR rate. 

DR rate coefficients are sometimes calculated using the Burgess (1965) general formula (GF) 
or versions of the GF as modified by Merts et al. (1976) and by Burgess & Tworkowski (1976). 
We have calculated the Burgess and Merts et al. rates using the formulae as given by Cowan 
(1981). We use oscillator strengths and excitation energies from Wiese, Smith, & Glennon (1966). 
Shown in Figure |5| for C V is the Merts et al. rate which is in fair agreement with our derived 
rate, though with decreasing Tg it goes to zero sooner than the experimental rate. Also shown 
in Figure ^ is the Burgess GF rate which is a factor of ~ 2.2 times larger than the experimental 
rate. The Burgess & Tworkowski (1976) modification of the GF yields a rate ~ 50% larger than 
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our derived rate. This is surprisingly good considering that their modification is meant for DR 
forming heliumhke, not hthiumhke, ions. Figure ^ shows for O VIII the GF rate which is ^ 38% 
larger than the experimental rate and, considering the expected accuracy of the Burgess formula, 
in fair agreement. The Burgess & Tworkowski rate is in good agreement. Also shown is the Merts 
et al. formula rate which is a factor of ~ 2.2 smaller than the experimental rate. Our C V and 
O VIII results strongly suggest that for AA^ = 1 DR it is not possible a priori to know which 
formula will yield a result closer to the true rate. 

For use in plasma modeling, we have fit the experimentally derived C V and O VIII DR rates 
using the simple formula (Arnaud & Raymond 1992) 

a(Te)=^-3/2^c,e-^»/*^^^^ (9) 

i 

Here Cj and Ei are, respectively, the strength and energy parameters for the ith fitting component. 
Best fit values are listed in Table |. For C V, the fit is good to better 1% for 2.6 x 10^ <Te< 10^ 
K. Below 2.6 x 10^ K, with decreasing Tg the fit goes to zero faster than the derived rate. But this 
error is insignificant for plasma modeling as the RR rate is ~ 3 orders of magnitude larger than the 
DR rate at these temperatures. For O VIII, the fit is good to better 1% for 7.1 x 10^ < Tg < 10* 
K. Below 7.1 X 10^ K, with decreasing Tg the fit goes to zero faster than the derived rate. This 
error is insignificant for plasma modeling as the RR rate is ~ 2 orders of magnitude larger than 
the DR rate at these temperatures. 



4. Summary 

We have presented a simple technique for obtaining DR rate coefficients from laboratory 
measurements of DR resonance strengths and energies. With this technique, we have derived DR 
rates for C V to C IV and O VIII to O VII using published resonance strengths and energies. 
Our derived rates are in good agreement with multiconfiguration, intermediate-coupling and 
multiconfiguration, fully-relativistic calculations as well as with most L5'-coupling calculations. 
Our rates are not in agreement with the recommended DR rates commonly used for astrophysical 
plasma modeling. We have used theoretical radiative recombination (RR) rates in conjunction 
with our derived DR rates to produce a total recombination rate for comparison with unified 
RR+DR calculations in LS coupling. Our results are not in agreement with undamped, unified 
calculations for C V but are in reasonable agreement with damped, unified calculations for O VIII. 
Also, neither the Burgess general formula, the Merts et al. formula, nor the Burgess & Tworkowski 
formula consistently yield a rate which is in agreement with our derived rates. This suggests that 
for An = 1 DR it is not possible to know a priori which formula will yield a rate closer to the true 
DR rate. We have also presented simple fitting formula of the experimentally derived DR rates for 
use in plasma modeling. 
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Table 1. Fit parameters for the experimentally derived C V to C IV and O VIII to O VII DR 
rate coefficient. The units for q are cm^ K^'^ and for Ei are eV. 





C V 




VIII 


i 


Ci 


Ei 


Ci Ei 


1 


3.02e-3 


246 


8.56e-3 497 


2 


1.75e-2 


296 


5.43e-2 632 
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Temperature (K) 

Fig. 1. — C V to C IV Maxwellian-averaged DR rate coefficients. The thin sohd curve is the 
integration of the experimental DR resonance strengths and energies from Kilgus et al. (1993) and 
Mannervik ct al. (1997). The short-dashed curve the recommended rate of Shull & van Steenberg 
(1982). The long-dashcd curve is the single-configuration, LS'-coupling results of Romanik 
(1988) and the long-dashed-dot curve of Bellantone & Hahn (1989). The triangles are the LS'- 
coupling calculations of Younger (1983); the circles the intermediate-coupling, multiconfiguration 
calculations of Badnell et al. (1990); and the squares the fully-relativistic, multiconfiguration 
calculations of Chen (1988). The thick solid line is the RR rate of Verner & Ferland (1996). 
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Fig. 2. — O VIII to O VII Maxwellian-averaged DR rate coefficients. Tlie tfiin solid curve is tfie 
integration of tfie experimental DR resonance strengths and energies from Kilgus et al. (1990). The 
short-dashed curve the recommended rate of Shull & van Stccnbcrg (1982). The long-dashed-dot 
curve is the single-configuration, LS-coupling results of Bcllantonc &: Hahn (1989) and the dotted 
curve the intermediate-coupling, multiconfiguration calculations of Pindzola et al. (1990). The 
thick solid line is the RR rate of Verner &; Ferland (1996). 
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Fig. 3. — C V to C IV Maxwellian-avcragcd RR+DR rate coefficients. The solid curve is sum of 
the RR rate from Verner k. Ferland (1996) plus the DR rate derived from the integration of the 
experimental DR resonance strengths and energies from Kilgus et al. (1993) and Mannervik et al. 
(1997). The dotted curve is the undamped, unified RR+DR calculations in LS'-coupling of Nahar 
& Pradhan (1997). 
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Fig. 4. O VIII to O VII Maxwcllian-avcraged RR+DR rate coefficients. The solid curve is sum 
of the RR rate from Verncr & Ferland (1996) phis the DR rate derived from the integration of the 
experimental DR resonance strengths and energies from Kilgus ct al. (1990). The dotted curve is 
the damped, unified RR+DR calculations in LS-coupling of Nahar (1999). 
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Fig. 5. — C V to C IV MaxwcUian-averaged DR total rate coefficients. The solid curve is the 
integration of the experimental DR resonance strengths and energies from Kilgus et al. (1993) and 
Mannervik et al. (1997). The dotted curve is the Burgess (1965) formula, the long-dashcd curve the 
Merts et al. (1976) formula, and the short-dashed curve the Burgess & Tworkowski (1976) formula. 
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Fig. 6. O VIII to O VII MaxwcUian-averaged DR total rate coefficients. The solid curve is the 
integration of the experimental DR resonance strengths and energies from Kilgus et al. (1990). The 
dotted curve is the Burgess (1965) formula, the long-dashed curve the Merts et al. (1976) formula, 
and the short-dashed curve the Burgess &: Tworkowski (1976) formula. 



